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Abstract — This study evaluated Earth Observing 1 (EO-1) 
Hyperion reflectance time series at established calibration sites 
to assess the instrument stability and suitability for monitoring 
vegetation functional parameters. Our analysis using three 
pseudo-invariant calibration sites in North America indicated that 
the reflectance time series are devoid of apparent spectral trends 
and their stability consistently is within 2.5-5 percent throughout 
most of the spectral range spanning the 12+ year data record. 
Using three vegetated sites instrumented with eddy covariance 
towers, the Hyperion reflectance time series were evaluated for 
their ability to determine important variables of ecosystem func- 
tion. A number of narrowband and derivative vegetation indices 
(VI) closely described the seasonal profiles in vegetation function 
and ecosystem carbon exchange (e.g., net and gross ecosystem 
productivity) in three very different ecosystems, including a hard- 
wood forest and tallgrass prairie in North America, and a Miombo 
woodland in Africa. Our results demonstrate the potential for 
scaling the carbon flux tower measurements to local and regional 
landscape levels. The Vis with stronger relationships to the C0 2 
parameters were derived using continuous reflectance spectra and 
included wavelengths associated with chlorophyll content and/or 
chlorophyll fluorescence. Since these indices cannot be calculated 
from broadband multispectral instrument data, the opportunity 
to exploit these spectrometer-based Vis in the future will depend 
on the launch of satellites such as EnMAP and HyspIRI. This 
study highlights the practical utility of space-borne spectrometers 
for characterization of the spectral stability and uniformity of 
the calibration sites in support of sensor cross-comparisons, 
and demonstrates the potential of narrowband Vis to track and 
spatially extend ecosystem functional status as well as carbon 
processes measured at flux towers. 

Index Terms — Earth Observing-1, eddy covariance, flux sites, 
forest, grassland, Hyperion, imaging spectroscopy, Pseudo-In- 
variant Calibration Sites, reflectance, remote sensing, seasonal 
carbon dynamics, spectral time series, vegetation function, wood- 
land. 
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I. INTRODUCTION 

W ITH the changing climate, it has become critical to 
understand land cover dynamics as ecosystems cycle 
through seasonal changes and respond to variable environ- 
mental conditions. The ability to monitor ecosystem carbon 
accumulation is of great interest when evaluating the carbon 
balance between various ecosystems. Long term satellite 
observations and land cover records are required for such 
environmental monitoring and change detection. However, 
multi-date satellite monitoring of the terrestrial environment 
requires adequate radiometric stability of the data to carry 
out biophysical and geophysical parameter surveys that are of 
sufficient sensitivity and accuracy, and are reproducible over 
time. 

The Hyperion 1 spectrometer carried on the Earth Observing 
One (EO-1) satellite has provided a rich high resolution data 
record over more than a decade (2001-2012+). Hyperion 
spectra cover the 400 to 2500 nm range with 242 overlapping 
spectral bands (196 of which are well calibrated) at approxi- 
mately 10 nm spectral resolution and 30 m spatial resolution, 
typically imaging a 7 km x > 100 km scenes. Use of these 
hyperspectral images to characterize terrestrial surface states 
and processes has rapidly increased since 2009, when the EO-1 
Hyperion archive was made available at no cost by United 
States Geological Survey (USGS) [1]. Hyperion data have 
several advantages over multispectral satellite systems: they 
provide information critical for atmospheric correction of top 
of atmosphere (TOA) radiances to derive surface reflectance, 
they enable the use of a broad array of spectral indices derived 
throughout the continuous spectra, and they provide the ability 
to simulate bands equivalent to broadband systems covering 
similar spectral ranges [2]. 

Hyperion data have been successfully used in discriminating 
among land cover types and species groups [3], deriving Vis 
[2], [4], estimating biophysical products such as forest canopy 
nitrogen [5], [6], primary production [6], forest canopy closure 
[7], vegetation fractional cover [8], [9], and canopy biophysical 
properties such as greenness, wetness and pigment content [10]. 

Earlier studies, describing the signal to noise ratio (SNR) of 
Hyperion’s radiance measurements in comparison to AVIRIS, 
the Landsat series and others [10]— [14], have found the data 
comparable to ~1990 Airborne Visible/Infrared Imaging Spec- 
trometer (SNR ~ 160:1 in the 0.4-1. 0 /im; and ~40:1 in the 

1 Any use of trade, firm, or product names is for descriptive purposes only and 
does not imply endorsement by the U.S. Government. 
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1.0-2. 5 /im). Aggregating the 10-nm Hyperion bands into 
simulated ETM+ multispectral bands has allowed radiometric 
cross calibration between the Hyperion and Landsat ETM+, 
or MODIS data, for comparison and conformation of the 
obtained results [15]. However, most of these studies analyzed 
individual images, or carefully selected scenes from the same 
season. In 2008, the EO-1 mission expanded the collection of 
time series for calibration and validation efforts, to assemble 
datasets essential for vegetation assessments. Currently there 
are collections of ten or more images at select sites with eddy 
covariance towers for examining the dynamics of various land 
covers at high spectral and spatial resolutions [1]. 

The goal of this study is to assess the potential of Hyperion re- 
flectance time series to describe seasonal changes in ecosystem 
functioning. Specific objectives include: 1) evaluation of the 
temporal stability of Hyperion’s reflectance measurements over 
three pseudo -invariant calibration sites (PICS) which have little 
seasonal variation in land cover, and 2) quantification of the 
relationships between Vis from Hyperion’s reflectance time 
series and ecosystem carbon accumulation (e.g. Net and Gross 
Ecosystem Production), measured over three validation flux 
sites (VFS) which have strong and very different seasonal leaf 
area index patterns. 

The incoming surface radiance signal is subject to several fac- 
tors that influence the stability and sensitivity of Hyperion (and 
any orbital) spectral data. These include sensor shortcomings 
such as relatively low SNR, and deficiencies in the informa- 
tion and algorithms utilized in initial calibration and processing 
streams (e.g., standards, measurements, models and accuracy 
of calibrations performed along the data processing chain). To 
study surface properties using imaging spectrometer data, ac- 
curate removal of atmospheric absorption and scattering effects 
(e.g., through implementation of atmospheric correction rou- 
tines) is necessary to derive surface reflectance from top of at- 
mosphere (TO A) radiances. These procedures account for the 
differences in solar illumination during the year at various lo- 
cations and correct for differences in atmospheric conditions at 
the time of acquisition [16]. However, they may also introduce 
variations into the reflectance data, due to differences between 
the modeled and real atmospheres at the time the data were col- 
lected [2], [16]. As a result, the cumulative error may lead to 
a variation in the reflectance signal of several percent in each 
spectral band, which could severely limit the sensitivity of the 
data, and/or introduce an apparent spectral trend. 

Vegetation properties are often assessed by computing Vis 
using a combination of two or more spectral bands. Vis are 
designed to relate more clearly to biophysical variables (e.g., 
chlorophyll, water content) than the original bands alone [17]. 
Also they tend to minimize variations not related to the variables 
of interest, such as illumination levels [4], [17]. Vis have been 
applied to the characterization of vegetation cover using broad- 
band systems, but have even more robust application using hy- 
perspectral data and specific narrowbands for the estimation of 
vegetation function (e.g., stress, water and light use efficiency, 
and chlorophyll fluorescence parameters) [5], [17], [19] and a 
range of leaf-level traits (e.g., pigments, cellulose and water 
content) [5], [1 8]— [21]. Roberts et al. [4] summarized Vis into 
three general categories of vegetation properties: structural, bio- 


chemical, and plant function/stress. Recent reviews of hyper- 
spectral Vis [4], [17] note that their use with reflectance time 
series remains to be thoroughly examined. 

II. Methods 

A. Study Sites 

To assess the stability of the Hyperion reflectance spectra we 
used three desert calibration sites and three vegetated validation 
sites. The vegetated sites were selected to provide a variety of 
vegetation types. 

1) Pseudo-Invariant Calibration Sites (PICS): The three 
desert Pseudo-Invariant Calibration Sites (PICS) are located 
in the western USA and include: Frenchman Flat (FMF), 
Ivanpah Playa (IP) and Railroad Valley Playa (RRVP). They 
are characterized by high reflectances, high spatial and temporal 
uniformity, high sun elevation, and minimal cloud cover. The 
PICS (Table I) are among the sites endorsed by the Committee 
on Earth Observing Satellites (CEOS, http://www.ceos.org/) 
as standard references for the post-launch calibration of 
space-based optical imaging sensors [22]. The FMF site is 
situated in the homogeneous section of the Frenchman Flat 
dry lakebed located NNE of Mercury, Nevada, USA on the 
Nevada Test Site range. It serves as The FED-based Spectral 
(FSpec) vicarious calibration test site [23]. The site is also 
included in the Aerosol Robotic Network (AERONET [24]), 
and is ideal for use with sensors with less than 300 m pixel 
size. The IP site is a dry-lake playa, located in the Prim Valley 
on the border of California and Nevada, USA and is large 
enough for use with mid spatial resolution sensors (~1 km). 
The RRVP site is located in a large dry lakebed in central 
Nevada, approximately 300 miles north of Las Vegas and 100 
miles east of Tonopah, and has a dry climate typical of the high 
desert of western USA [25]. Because of its large size, RRVP 
is suitable for sensors with even larger footprints (1-10 km), 
and is automated with instrumentation used extensively for the 
vicarious calibration of terrestrial imaging sensors covering 
the visible (VIS), near-infrared (NIR) and shortwave infrared 
(SWIR) wavelength ranges [26], [27]. 

The surface layers and composition of all three sites are rel- 
atively smooth and spatially homogeneous, consisting of com- 
pacted clay-rich lacustrine deposits [22]. All sites suffer from 
the presence of iron (Fe 3+ ) absorption features, which effects 
the spectral properties and is a typical characteristic of playas in 
this region of the USA [22], [25]. Compared to RRVP, the sur- 
face layers at FMF and IP contain less loose sand and evaporate 
salts, and are therefore more stable and spatially uniform [22], 
[25] — [28]. 

2) Vegetated Validation Flux Sites (VFS): The selected Vali- 
dation Flux Sites (VFS) are located in Konza Prairie near Man- 
hattan, Kansas, USA, the Duke Forest near Durham, North Car- 
olina, USA, and Mongu, Zambia in Africa. Their characteristics 
are also described in Table I. The vegetated sites were part of 
the FLUXNET network, and had instrumented towers providing 
meteorological and eddy covariance data (Table I). The Konza 
and Mongu VFS are a part of the network of Earth Observing 
System (EOS) Land Validation Core Sites for land product val- 
idation. 
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TABLE I 

Description of Calibration and Validation Study Sites: Location and Land Cover Type/Species Composition 


Site Name 

Center 

Latitude 

Center 

Longitude 

Elevation 
(in asf) 

Land Cover and Vegetation Type 

Frenchman 
Flat, NV, USA 

36.80928 

-1 15.93479 

940 

Play a, no vegetation 

Ivanpah Plava 
CA, USA 

35.5692 

-1 15.3976 

813 

Play a, no vegetation 

Railroad Valley 
NV, USA 

38.50 

-1 15.69 

1435 

Play a, no vegetation 

Duke Forest 
NC, USA 

Konza T KZ 
USA 

Mongu : Zambia 

35.973583 

39.080002 

-15.4377 

-79.100433 

-96.559998 

23.2527 

170 

320-440 

1080 

Mixed Hardwoods: Carya tomentosa, C glabra , 
Liriodendron tulipifera, Liquidambar styraciflua, 
and oaks (Q. alba> 0 michauxii, Q. phellos }; 
Loblolly plantation: Pi mis laeda L 
Tallgrass prairie. C4 grasses: Andropogon 
gerardii, Sorghastrun nutans, Panicum virgatum, 
Schizachyrium scopari nm 
Miombo Woodland; Brachystegia bake ri ana 
Brachystegia spiciformis (95/5 %) 


t asl - above sea level 


The Duke Forest (North Carolina, USA) study site includes 
two adjacent ecosystems, with flux towers operating during 
2001-2010 by Duke University. One of the towers was located 
in an even-aged 20+ year-old loblolly pine (LP) plantation and 
the other in a mixed hardwood stand (HW). The two forests 
differ in predominant vegetation cover type and canopy struc- 
tures, while experiencing nearly identical climatic and edaphic 
conditions. The LP site was established in 1983 following a 
clear cut and a burn events [28], [29], and is comprised pri- 
marily of P. taeda with some emergent Liquidambar styraciflua 
L. and a diverse sub-canopy. The mixed HW mature forest is 
a 90- 11 0-year old stand, dominated by oaks ( Quercus alba 
L., Q. michauxii Nutt., Q. phellos L.) and hickories ( Carya 
tomentosa (Poir.) Nutt., C. glabra (P. Mill). Sweet.), with 
some yellow poplar ( Liriodendron tulipifera L.) and sweetgum 
( Liquidambar styraciflua L.) [28], [29]. 

The Konza Prairie site (Kansas, USA, Table I) includes two 
flux towers located in grasslands managed under one- and four- 
year burn cycles [30]-[32]. The vegetation is primarily (>90%) 
native tallgrass (TG) prairie species dominated by perennial 
C 4 grasses such as Andropogon gerardii , Sorghastrun nutans , 
Panicum virgatum and Schizachyrium scoparium. However, nu- 
merous sub-dominant grasses, forbs and woody species con- 
tribute to the high floristic diversity of the Konza Prairie site 
[30], [31]. It has a continental climate characterized by warm, 
wet summers and dry, cold winters. The mean annual precipi- 
tation of 835 mm is sufficient to support woodland or savanna 
vegetation. Consequently, drought, fire and grazing are impor- 
tant in maintaining this topographically complex grassland [31], 

[32] . 

The African validation site near Mongu, Zambia (Table I) is 
monitored by one eddy covariance tower located in a forested 
area [33], [34] where flux measurements have been acquired 
since 2008. The site is dominated by hardwood trees and re- 
ferred to as Miombo woodland (MW) with a fractional cover of 
65% as measured by Scholes et al. [35]. Grass cover is minimal 

[33] with the remaining surface area consisting of bare aeolian 
sand. Canopy cover heterogeneity is largely the result of distur- 
bances associated with subsistence forestry, grazing, and fire, 
which contribute to the patchiness of the land surface. Based 
on a 20-yr average (1973-93) mean annual rainfall for Mongu 


was 879 mm, with 94% of this amount received during the wet 
season months of October-March [34]. 

B. Data Collections and Processing 

1) EO-1 Hyperion Time Series Collection: The Hyperion 
image collection in the 2001-2011 timeframe over the study 
sites is comprised of more than 350 images with varying extents 
of cloud cover. The comparisons and analyses were restricted 
to cloud- free (< 20 %) data with a maximum of five degrees 
difference in sensor viewing geometry. The optimal Hyperion 
images (Table II) were selected and downloaded as calibrated 
at-sensor radiances from the USGS (EarthExplorer, http://earth- 
explorer.usgs.gov/), and processed using ENVI [36], ERDAS 
Imagine [37] and PRISM [38]. 

The TOA radiances measured by Hyperion were converted 
to surface reflectances, using the commercially-available soft- 
ware, Atmosphere CORrection Now (ACORN) [39], designed 
for pushbroom imaging spectrometers with cross-track spectral 
calibration variation. ACORN uses the MODTRAN4 [40] 
radiative transfer model, constrained by the elevation and the 
observation geometry, to explicitly simulate the absorption and 
scattering effects of atmospheric gases and aerosols, in order 
to produce apparent surface reflectance. A key characteristic 
of ACORN is the use of full spectral fitting methods to solve 
for the overlap of absorptions between water vapor and liquid 
water in surface vegetation [39]. From calibrated Hyperion 
radiances, the ACORN routine (mode 1 pb) produced a high 
spectral resolution apparent surface reflectance image cube and 
a water vapor single-band image for each acquisition (Table II) 
estimated on a pixel-by-pixel basis using the absorption bands 
at 940 and 1150 nm. Both transmitted radiance and atmo- 
spheric reflectance were calculated for each pixel using the 
derived water vapor, atmospheric pressure, site elevation, and 
aerosol optical depth estimates. Apparent surface reflectance 
was derived using either mid-latitude summer or winter atmo- 
spheric models [39] depending on the acquisition date and site 
location. ACORN offers artifact (e.g., spectrally incoherent 
noise) suppression options, which were not employed in order 
to preserve original spectral variability. Also, the images were 
not geographically or geometrically rectified. 
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TABLE II 

EO-1 Hyperion Images Used in This Study 


Site Name 

Cloud- free 
Images 

[number (year)l 

Analyzed Scenes by Acquisition Date 
[doy/year] 

Frenchman Flat, USA 
Ivanpah PI ay a, USA 
Railroad Valley, USA 

6 (2003-2009) 
6 (2002-2005) 
14 (2001-2008) 

12/2009. 178/2003. 199/2008. 219/2007, 222/2008, 270/2007, 
337/2008. 

68/2003.81/2002, 148/2003, 169/2005, 180/2003.244/2003. 
12/2006,41/2006, 120/2007, 135/2007, 165/2001, 174/2004, 
176/2008. 181/2001, 197/2001,209/2007, 234/2006, 
272/2005, 284/2005, 288/2007. 

Duke (HW and Pines), USA 

7 (2008-2009) 

6, 34, 162, 180, 203, 290, 300. 

Konza Prairie, USA 

7 (2008-2009) 

13, 28, 202, 205, 26 ! t 266, 302 

Mongu, Zambia 

25 (2008-2009) 

22,40. 107, 115, 141, 145, 158, 163, 171, 173, 176, 181. 191. 
194. 199, 201, 204, 225, 230,253. 284, 291, 296, 302, 365. 


TABLE III 

Subset of Hyperion 224 Bands to 171 at PICS and 159 VFS 


Array 

Pseudo Invariant Cal. Sites (PICS) 

Vegetative Sites 

Bands 

Wavelengths (nm) 

Bands 

Wavelengths (nm) 

VN1R 

8-57 (49) 

426-926 

16-57(49) 

498-926 

SW1R1 

79-120(41) 

993;-l 346 

78-120(41) 

993:-I346 

SWIR2 

129-165 (36) 

1518-1800 

137-164(27) 

1518-1790 

SW1R3 

179-224(45) 

1942-2395 

198-213(42) 

2133-228 5 

J otal 

LULL 


(159) 



Band subsetting was used to remove uncalibrated and over- 
lapping bands, and bands with low SNR and high variability 
(e.g., bands adjacent to water absorption features), resulting in 
subsets of 171 bands for PICS and 159 for VFS, covering por- 
tions of the VIS, NIR and SWIR regions (Table III). The smaller 
spectral band subset at the VFSs is due to the stronger effect of 
canopy water absorption features on the adjacent bands, causing 
errors in the atmospheric correction procedure. Bands in wave- 
length ranges around 1430 nm and 1930 nm were also elimi- 
nated due to strong disturbance by atmospheric water vapor. 

Local spatial statistics were used to discern spatially homo- 
geneous areas for spectral extraction at the PICS. We evaluated 
spatial uniformity for all dates by calculating the Getis Ord (Gi) 
statistics and Moran’s I index in ENVI [41]. A cluster of pixels 
with high (or low) digital count values is indicated by largely 
positive (or negative) Gi values [36], [41]. Moran’s I index, 
a measure of local spatial uniformity, ranges between + 1 and 
— 1, where + 1 defines spatially correlated pixels (i.e. those with 
clustered data values) and — 1 spatially uncorrelated data (i.e., 
pixel values with no clustering). We also calculated mean re- 
flectance at each band and used the coefficient of variation (CV) 
and the ratio of the standard deviation (SD) over the mean (x) 
[22], to characterize the radiometric and temporal stability for 
each band. 

We identified the position of the FLUXNET sites on the im- 
ages using available maps, photos and images. Spectra were ex- 
tracted, from corresponding areas, from 20-30 pixels within the 
Hyperion images. We calculated mean reflectance at each band 
and used the coefficient of variation (CV) and the ratio of the 
standard deviation (SD) over the mean (x), to characterize the 
temporal variation for each band. The slight shifts in the overlay 
of the pixels of spectral extraction on the ground unavoidably 


introduce statistical uncertainty and random variation, which 
maybe more pronounced at the Mongu site, where the wood- 
land canopy cover is more sparse and variable. We chose to re- 
tain the nearly full-spectral range of Hyperion data, in order to 
test the ability of the band depth analysis and spectral Vis to de- 
tect changes in vegetation function and to track C0 2 seasonal 
dynamics measured at the FLUX towers. 

2) Eddy Covariance Data at Validation Flux Sites (VFS): 
Net Ecosystem Production (NEP) provides a comprehensive 
measure of ecosystem net carbon accumulation which is defined 
as the carbon accumulation within ecosystems and is the differ- 
ence between gross primary production (GPP) and ecosystem 
respiration (Re) [42]. NEP is measured with eddy covariance 
methods using sonic anemometers mounted on towers above the 
canopies [42]. For each VFS, we assembled a set of data (e.g., 
NEP and associated meteorological parameters such as precipi- 
tation and temperature) covering the timeframe of the Hyperion 
image acquisitions. For the analysis, the flux data and the Hy- 
perion reflectance were paired by date of acquisition. Data for 
the Mongu site were downloaded from the CarboAfrica web site 
(http://www.carboafrica.net/data_en.asp, PI for Mongu Dr. W. 
Kutsch), and flux data from the other sites were provided by 
their tower managers Dr. Brunsell (KNZ-LTER) for the Konza 
Prairie site, and Dr. Kimberly Novick (The School of the Envi- 
ronment at Duke University) for the Duke Forest sites. 

Midday flux averages were calculated as the average of the 
observations collected between 11:00 AM and 1:00 PM local 
time. The data for Mongu and Konza were not gap-filled. The 
Duke NC data were gap-filled using methods described in Stoy 
et al. 2006 [43]. 

NEP was partitioned into Gross Ecosystem Production (GEP) 
and Ecosystem Respiration (Re). For the Duke and Konza sites, 
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daytime Re was calculated from relationships developed be- 
tween nighttime NEP and air temperature [44], [45]. For the 
Mongu site, air temperature alone did not provide a good de- 
scription of variation in Re, so the Re model used both air tem- 
perature and soil moisture. GEP was computed as the sum of 
NEP and Re. 

3) Reflectance Time Series for Validation Sites: The Hy- 
perion time series included multiple scenes in order to con- 
struct a phenology season, with many images collected over 
multiple years. Due to the non-continuous temporal resolution 
(i.e. tasked and targeted image acquisition) of Hyperion datasets 
used in this study (Table II), annual documentation of seasonally 
representative vegetation phenological changes was not pos- 
sible. However, both spectra and the corresponding eddy co- 
variance CO 2 flux measurements were organized in seasonally 
successive chronological order based on the day of year (DOY), 
regardless of acquisition year. 

4) Continuum Removal and Vis From High Spectral Resolu- 
tion Spectra: We applied the ‘band depth analysis of absorption 
features’ method [46] that enhances and standardizes known 
spectral chemical absorption features. Continuum removal nor- 
malizes reflectance spectra to allow comparison of individual 
absorption features from a common baseline [38], [46], [47]. 
The continuum is a straight line segment fitted over the hull 
of the spectrum, thus connecting the local maxima [46]. The 
continuum-removed reflectance is obtained by dividing the re- 
flectance value in the absorption minimum by the reflectance 
level of the continuum line at the corresponding wavelength. 
The Processing Routines in IDL for Spectroscopic Measure- 
ments (PRISM) provide a uniform and automated approach for 
obtaining continuum removed (CR) Hyperion reflectances, with 
spectral profile output having values between 0 and 1 for en- 
hancing the selected features [38], [46]. 

Six known chemical absorption features have been associ- 
ated with foliar chemistry, including the VIS chlorophyll ab- 
sorption features (470-520 nm, 550-750 nm), and SWIR ab- 
sorption features (1110-1285 nm, 1630-1790 nm, 2005-2195 
nm, and 2220-2380 nm) detecting water, lignin and cellulose 
[38], [46], [48]. We applied the band depth analysis to study 
the seasonal changes in feature depth (FD) in the well-defined 
chlorophyll absorption feature (between 500-800 nm). Feature 
depth was computed as 


band centers and 10 nm resolution, we generated 76 high res- 
olution Vis: 18 for structure (e.g., LAI and biomass), 31 for 
biochemical content (e. g., chlorophyll, 18; carotenoids, 2; an- 
thocyanins, 2; photo synthetic pigments, 5; water, 3; and lignin 
and cellulose, 1); 27 for plant function (e.g., stress, 12; light use 
efficiency, 4; and chlorophyll fluorescence, 11). We were not 
able to generate Vis that required the use of blue wavelengths 
below 450 nm (e.g., SIPI for pigment content), or wavelengths 
between 1350-1500 nm (for canopy structure, water content) 
and 1950-2100 nm (for lignin and cellulose content), due to 
the lower data quality at these wavelengths. These 76 useful 
Vis included band ratios of derivative (typically first derivative) 
spectra, which have been demonstrated to be strongly related to 
chlorophyll concentration, chlorophyll fluorescence and forest 
stress [19], [20]. The first derivatives (D) of reflectance were 
calculated using the formula: 


(B-Mi+1) - R A(j-i)) 

(A (i+ i) - A (i _ 1} ) 


( 2 ) 


where D\[ is the first derivative and R\[ is the mean reflectance 
value at Ai wavelength. The derivative transformation provides 
spectral data that are less sensitive to the effects of variable ir- 
radiance and the influences of canopy structural and bidirec- 
tional properties. To compare the above high spectral resolution 
Vis with broadband Vis we simulated the bands of Landsat TM 
(lower spectral resolution), and calculated the most widely used 
Normalized Difference Vegetation Index (NDVI) and the En- 
hanced Vegetation Index (EVI). 

5) Statistical Analysis : Analysis of Variance routines and de- 
scriptive statistics [49] are used to assess temporal spectral uni- 
formity and/or Vis at each site. The resulting mean (x) and stan- 
dard deviation (SD) statistics are reported as plots for each site 
and summarized for comparison between the sites. To examine 
the temporal variability in the spectral properties of VFS, the 
means and standard deviations for each of the sites were plotted 
and discussed further in the results section. Also at the VFS, 
analyses of variances and regression procedures were used to 
relate the Hyperion Vis to NEP and GEP tower measurements 
[49]. Spectra were analyzed individually for each VFS, as well 
as collectively as a multi-site dataset. 


FD = (1 - CR), (1) 


III. Results 


where CR is the continuum removed reflectance at 675 nm, so 
that increasing amounts of chlorophyll absorber yield increasing 
FD. We selected this strong feature because it is in a high Hy- 
perion S/N region, away from atmospheric water absorption 
bands, and thus less subject to noise or residual atmospheric 
contamination. 

Recent review articles [4], [17] provide a summary of the 
utility of Vis derived from high resolution (< 10 nm bands) re- 
flectance and derivative spectra [19], [20]. At the VFS, we tested 
the ability of a subset of reflectance Vis applicable for agricul- 
ture and forestry to track the dynamics in vegetation C0 2 as- 
similation parameters (e.g., NEP and GPP). Using the full spec- 
tral range of Hyperion, and adjusting the formulas to Hyperion’s 


A. Seasonal Dynamics of Hyperion Reflectance at 
Pseudo-Invariant Calibration Sites (PICS) 

The Gi statistics and Moran’s I index describe the spatial 
uniformity of the PICS (Figs. 1 and 2, Table IV). The FMF 
site offers an approximately 2x2 km area with relatively uni- 
form high reflectance distribution, as determined by the Gi sta- 
tistics (Figs. 1(a), 2(a)). The homogenous area at IP was ap- 
proximately 2x4 km (Figs. 1(b), 2(b)). RRVP had an area of 
approximately 6x10 km with relatively uniform reflectance 
(Figs. 1(c), 2(c)). The highest reflectance values were measured 
across the 650-2000 nm region at FMF, followed by IP and 
RRVP (Figs. 1, 2). At RRVP, we observed that certain areas of 
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Fig. 1. Pseudo invariant calibration sites (PICS): (a) Frenchman Flat (FMF), (b) Ivanpah Playa (IP), and (c) Railroad Valley (RRVP). The center of the region 
(30^10 pixels) of reflectance extraction is marked with a cross. Two portrayals are shown for each site from a summer acquisition: natural color composites (fig. 
(a)l, (b)l, (c)l; RGB: 651,549,447 nm), and the Getis Gi statistics calculated to indicate homogeneous pixel clusters (fig. (a)2, (b)2, (c)2). 
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Fig. 2. Temporal variation (X d= SD) in the Hyperion reflectance at select wavelengths: (a) Frenchman Flat (FMF, 7 images acquired in 2005-2008), (b) Ivanpah 
(IP, 6 images acquired in 2003-2005) and (c) Railroad Valley (RRVP, 15 images acquired in 2001-2008). The seven Hyperion bands selected to represent the 
seasonal reflectance trends are centered at: 447, 549, 651, 854, 1003, 1679 and 2204 nm. 


the site are brighter than others. Since the site is completely de- 
void of vegetation, this contrast can be explained by soil mois- 
ture effects and variations in the mineral composition of the 
playa’s surface. 

The three calibration sites (Fig. 3) can be characterized as 
having high mean reflectance values, across a wide spectral 


range, which agrees with previous studies conducted over the 
same sites using field radiometers and/or space-borne multi- 
spectral sensors [12], [14], [37]. Mean reflectance in the blue re- 
gion (450-550 nm) was highest at FMF (25-40%, ±2. 3-2. 5%, 
Fig. 3(a)), lowest at IP (8-12, ±2. 1-2.8%, Fig. 3(b)), and inter- 
mediate at RRVP (18-28, ±2-2.5, Fig. 3(c)). 
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TABLE IV 

Hyperion Reflectance and Moran I Statistics for Seven Wavelengths at These Pseudo Invariant Calibration Sites (PICS) 


Wavelength 

(nm) 

Frenchman Flat 

Ivan pah Playa 

Railroad Valley Playa 

Reflectance (%) 

Moran's J 

Reflectance (%) 

Moran's 1 

Reflectance (%) 

Moran's 1 

Mean 

SD 

CV 

(±d* 

Mean 

SD 

CV 

(±i)* 

Mean 

SD 

CV 

(± 1)* 

447 

0.21 

0.020 

0.095 

0.97 

0.17 

0.021 

0.122 

0.94 

0.20 

0.019 

0.095 

0.86 

549 

0.38 

0.024 

0.063 

0.98 

0.29 

0.031 

0.105 

0.97 

0.34 

0,025 

0.074 

0.87 

651 

0.50 

0.028 

0.056 

0.97 

0.38 

0,027 

0.071 

0.97 

0.37 

0.027 

0.073 

0.86 

854 

054 

0,032 

0.059 

0,98 

0,44 

0.025 

0,056 

0.96 

0,39 

0.028 

0.072 

0,95 

1003 

0.56 

0.037 

0.066 

093 

0.45 

0024 

0.053 

0.88 

0.40 

0,026 

0.065 

0.83 

1679 

0.58 

0.038 

0.066 

0.90 

0.5 

0.013 

0.026 

0.87 

0.42 

0.033 

0.079 

0.81 

2204 

0.48 

0.039 

0.081 

0.91 

0.42 

0.030 

0.071 

0.89 

0.32 

0.034 

0.106 

0,88 


+ Mean (x), Standard Deviation (SD), Coefficient of Variation (CV = SD/x) 

* Moran’s I index expresses local spatial uniformity, ranges between +1 and -1 (+1 = strong positive spatial autocorrelation, 
0 = spatially uncorrelated data, — 1 = strong negative spatial autocorrelation) 


However, in the green and red region (550-700 nm), reflec- 
tivity was highest at FMF (40-45, ±2. 5-2. 8%), intermediate at 
IP (35^0, ±2.4-2. 8%), and lowest at RRVP (30-40, ±2.5%). 
The same trend continued in the NIR (700-2000 nm): FMF 
50-55% (±3. 1-7.2%), IP 42^18% (±1.5-5%), and RRVP 
37-40% (±2.3— 5.8%); and in the SWIR (2000-2450 nm): 
FMF 48-52% (±3-18%), IP 40-48% (±2-6.5%), and RRVP 
32-45% (±3-7.5%). The highest SDs were associated with 
artifacts due to imperfections in the atmospheric correction 
procedures for scattering in the blue region and for absorption 
by water vapor near absorption features at 930 nm, 1450 nm 
and 1950 nm. In the 2360-2450 nm region, where the derived 
reflectance is less stable, due to lowered SNR of the instrument, 
the variation in the moisture content in the playa caused greater 
overall impact and the SD were consistently higher at all sites. 

The temporal variation in the reflectance characteristics of 
the three PICS are represented by individual bands in the blue, 
green, red, NIR, and SWIR regions (means and SD, Fig. 2). 
While the magnitude of the reflectance differed per spectral 
band, all reflective regions exhibited relatively uniform and 
flat temporal profiles, devoid of consistent temporal trends 
(Fig. 2). Lower reflectance values (with higher SD) across all 
spectral bands occurred during winter months at FMF (DOY 
12, 337) and RRVP (DOY 12, 288), which could be attributed 
to lower incoming radiation, higher shadow component, local 
illumination effects not accounted for in the atmospheric cor- 
rection model, and/or increased soil moisture at these playas. 
Fig. 3(a), (b), (c) and Table IV summarize the long-term spectral 
reflectance responses at the three PICS. The statistics (x, SD, 
CV) summarize band responses across time (Fig. 3, Table IV). 

The reflectance time series at the PICS were devoid of ap- 
parent spectral trends, allowing us to have sufficient confidence 
in the data to describe the temporal dynamics at the VFS sites. 

B. Seasonal Dynamics in Hyperion Reflectance Parameters 
and NEP at Validation FLUXNET Sites (VFS) 

Seasonal changes in Hyperion reflectance spectra over the 
forested Duke LP (Fig. 4(a)) and HW sites (Fig. 4(b)), the 
Konza prairie TG site (Fig. 4(c)), and the Mongu MW site 
(Fig. 4(d)) displayed typical vegetation reflectance trends: 


relatively low VIS reflectance due to high absorption by leaf 
photosynthetic pigments, with a peak ~550 nm, relatively 
bright NIR reflectance (700-1300 nm) because absorption is 
weak but scattering is strongly expressed by canopy structural 
components and the contrast between indices of refraction for 
water-rich cells and intercellular air spaces, and strong absorp- 
tion in the SWIR (1450-2450 nm) by leaf water, cellulose, 
lignin, and other plant materials causing vegetation to appear 
relatively dark. 

Hyperion’s reflectance time series captured seasonal dif- 
ferences among the site-specific VFS carbon dynamics. Red 
reflectance wavelengths (~700 nm) had the largest variation 
(large CV) during the year at both the grassland and deciduous 
hardwood sites where CV > 30% (HW, 50%; TG, 47%; and 
MW, 31%), but CV was less variable at the pine stand (LP, 
26%). The largest spectral variations were observed in the 
blue and red VIS spectrum (Table V). For blue wavelengths 
(450-500 nm), there were no statistically significant differences 
among the seasonal spectra for any VFS, which is attributed 
to the very large CV (70%-128%) in this region, due to the 
lower S/N of Hyperion’s first detector and imperfections in 
the correction for atmospheric scattering. Due to a pronounced 
wet/dry seasonality, the time series of red wavelengths at 
Mongu (Fig. 4) best tracked fluxes throughout the growing 
season. 

Fig. 5 depicts the temporal profile of the continuum removal 
analyses for reflectances in the chlorophyll absorption features 
at each VFS. The normalized reflectances have values between 
0 and 1, enhancing the feature and allowing feature depths (FD, 
measured at 675 nm) to be computed for temporal comparisons 
at each site, and among sites. At the Duke LP site, the FD did not 
vary significantly (during the year (mean FD = 0.81-0.84, p > 
0.25), while at the Duke HW site. FD dramatically changed from 
0.23 (DOY 6) to 0.90 (DOY162). The FD also varied greatly 
throughout the year at the Konza TG site, from 0.81 (DOY 205) 
to 0.05 (DOY13) and at the Mongu MW from 0.83 (DOY 22) to 
0.33 (DOY 253). At the Konza TG prairie site, the very low FDs 
(0.05-0.16) can be explained by the dormancy and complete 
dieback of the above ground tallgrass biomass. 

Attempts to analyze the dynamics in the SWIR absorption 
features (2220-2380 nm) associated with lignin, cellulose and 


CAMPBELL et al.\ EO-1 HYPERION REFLECTANCE TIME SERIES AT CALIBRATION AND VALIDATION SITES 


283 


(a) 


(b) 


= 

fZ 

E 


— reflectance 

- - - standard deviation 


0,75 




Z H r ' bv 0 . 5-1 

i / , r*-" 

| ,( W 


* 0.25 
■ 0,2 

■0.15 

-0J 

■0,05 

■0 


-Teflictance 
- - - standard deviation 


0.75- 




r' /V.. 

• U 


0.25 
| 0.2 
rO.I5 

I 

0.05 

0 


S o a o o o o o 

*ri O no O iri o to 

r . C '^ Tr - C ' M ^ 

- - - - rt (S 

Wavelength (nm) 


(c) 

l 

0 75 
0.5 


-reflectance 


0,25 




- standard deviation 


0-3 8 


-r^rsN 

^..AvP 


■ 0.15£ 

-0.1 <1 

T5 

b . 

0,05 -g 

a 

0 ^ 


<rt c trt c o 

On SM f* O* c-a Tt 

^ - - - n W 


Fig. 3. Mean reflectance (%, solid line) and standard deviation (%, dashed line) for time series for: (a) Frenchman Flat (FMF, n = 7), (b) Ivanpah Playa (IP, 
n = 6), and (c) Railroad Valley Playa (RRVP, n = 15). Means and standard deviations are listed in Table IV. 
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Fig. 4. Temporal variation in the Hyperion reflectance acquired over the validation sites: (a) Duke (loblolly pine, LP), (b) Duke (hardwoods, HW), (c) Konza 
(tallgrass, TG), and (d) Mongu (Miombo woodland, MW). Each line represents the mean site reflectance from a single day, designated in the labels as “DOY year”. 


nitrogen contents were not successful, due to low SNR which 
obscured the features and precluded the band depth analysis. 

Hyperion narrowband Vis were calculated using contin- 
uous reflectance spectra or reflectance derivative (PRI1, PRI4, 
Dmax), and included wavelengths associated with chlorophyll 
content (e.g., PRI4, G94 and G32) and/or chlorophyll fluo- 
rescence (e.g. Dmax/D704, Dmax/D( max+12 )) (Tables VI, 
VII). Using a regression approach, we established that many 
of the Vis were closely related to NEP and GEP. The Vis 
with the strongest relationships to NEP and GEP are listed in 
Table VI, along with the simulated Landsat broadband Vis 
(e.g., NDVI, EVI) for comparison. The set of best performing 
Vis for estimating NEP varied among the sites. At the Duke LP 
site, the most successful narrowband VI for association with 
NEP was the DPI (r 2 = 0.91) a derivative index using three 
wavelengths, followed by linear relationships obtained with 
two variations of the photochemical reflectance index: PRI4 


(r 2 = 0.85, Fig. 6(a)) and PRI1 (r 2 = 0.84). At the Duke HW 
site, the performance of PRI4 was non-linear but performed 
best (r 2 = 0.94, Fig. 6(a)), followed by Dmax (r 2 = 0.83). 
The PRI4 performed well at both LP and HW (Fig. 6(a)), how- 
ever it saturated at values >0.2, not being able to effectively 
discriminate for NEP > 20 ^mol m -2 s -1 . 

At the TG Konza site, however, the most successful 
Vis for relating best to NEP were the narrowband G94 
(r 2 = 0.95, Fig. 6(b), Table VI) and Dmax (r 2 = 0.93). 
G94 uses narrowband wavelengths in the NIR and red 
(r 2 = 0.95, Fig. 6(b)) and it saturated at G94 > 0.7 for 
NEP > 27 fimol m -2 s -1 . At the Mongu MW site, the most 
successful Vis for relating to NEP were G32 (r 2 = 0.83, 
Fig. 6(c)), Dmax/D704 (r 2 = 0.79), and Dmax (r 2 = 0.77). 
G32 combines wavelengths from the blue, red and NIR 
regions (Fig. 6(c)). At values of G32 <0.3 this parameter 
saturated for NEP < 5 ^mol m -2 s -1 . 
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TABLE V 

Seasonal Values lor ol Hyperion Reflectance at Validation FLUXNET Sites (VFS) 


Wavelength 

(nm) 

Duke, NC, USA 

Konza, KS, USA 

Mongu, Zambia 

Loblolly pine (LF) 

Hardwoods (HW) 

Tall Grassland (TG) 

Miombo Woodland (MW) 

Mean 

St. Dev 

CV 

Mean 

St, Dev 

CV 

Mean 

St. Dev 

CV 

Mean 

St. Dev 

CV 

447 

0.01 8m 

0.019 

1. 11 

0 01 Im 

0.008 

0.72 

0.035m 

0.044 

1.27 

0,013m 

0.017 

1.28 

549 

0,035f 

0,010 

0.29 

0.041 

0.010 

0,25 

0,060 

0,024 

0,40 

0.054 

0.015 

0,27 

651 

0.026 

0.007 

0.26 

0.039 

0.019 

Q.50 

0.073 

0.035 

0.48 

0.060 

0.019 

0,32 

854 

0.215 

0,026 

0,12 

0.259 

0.076 

0,29 

0,275 

0,076 

0,28 

0.238 

0 030 

0,13 

1104 

0.189 

0,031 

0,16 

0.232 

0.038 

0,16 

0,253 

0,033 

0,13 

0.253 

0.034 

0,13 

1679 

0.086 

0.01 1 

0.13 

0.155 

0.048 

0.31 

0.238 

0.037 

0.16 

0.206 

0.029 

0.14 

2204 

0.025 

0,006 

0.26 

0.070 

0.041 

0,59 1 

0.128 

0,041 

0.32 

0.108 

0.023 

0,21 


ns — no statistically significant differences among means by day of year (based on ANOVA) 

t when not noted as ns, there were statistically significant differences among means by day of year at p < 0.05 level 



Wavelength (ran) Wavelength (nm) 

Fig. 5 . Temporal variation in the mean normalized reflectance (continuum removed) from the chlorophyll absorption feature, (550-750 nm) illustrating the spectral 
dynamics at each validation site by DOY: (a) Duke LP (n = 6 days), (b) Duke HW (n = 6 days), (c) Konza TG (n = 6 days, each line is an average of spectra 
from areas under 1 and 4 year bum), and (d) Mongu MW (n = 15 days). 


At all VFSs, the Vis best related to NEP were also strongly 
related to GEP, with strongest relationship for the Konza TG site 
(Table VI). The broadband EVI also performed well at the Duke 
HW site (NEP, r 2 = 0.84), Konza TG site (NEP, r 2 = 0.83; 
GEP, r 2 = 0.89) and at the Mongu MW site (NEP, r 2 = 0.73; 
GEP,r 2 = 0.88). Fig. 7 depicts the best performing Vis for NEP 
across all sites. Most of these relationships were curvilinear 
(Fig. 6, Table VI). Hyperion narrowband Vis with the stronger 
relationship to CO 2 parameters (Tables VI and VII) were calcu- 
lated using continuous reflectance spectra or reflectance deriva- 
tive (e.g. PRI1, PRI4, Dmax), and included wavelengths associ- 


ated with chlorophyll content (e.g., PRI4, G94 and G32) and/or 
chlorophyll fluorescence (e.g. DPI, Dmax/D704). 

To assess the potential of using a common high resolution 
VI for tracking CO 2 seasonal dynamics across multiple sites 
and ecosystems, we combined the data from all VFS (n = 52, 
Table VII, and Fig. 8). The derivative parameters Dmax (as- 
sociated with vegetation stress) and DP22 (associated with 
chlorophyll fluorescence) exhibited a strong relationships to 
NEP and GEP (r 2 = 0.65-0.73 and r 2 = 0.74-77, respec- 
tively), as did the normalized difference water index (NDWI, 
r 2 = 0.74 to NEP and r 2 = 0.67 to GEP, Table VII). These 
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TABLE VI 

Hyperion Vegetation Indices (Vis) Having Strongest Relationships (r 2 ) to C0 2 Uptake Parameters (NEP and GEP), by Study Site. Hyperion’s 
Best Performing High Spectral Resolution Vis are Listed First. The Broadband Landsat Parameters are Listed for Comparison in Italics 


Index 

Bands (nm) 

to NEP 

r 2 to GE1’ 

L/NL* 

a. Loblollv Pine (LPL Duke. NC (n=7) 

DPI 

(D683* D7 L 0)/(D697* D697) 

0.91; p<0.01 

0.44, ns 

NL 

PRI1 

(R53 I-R570)/(R53 1-R570) 

0.84; pO.Ol 

0.73; p<0.0l 

L 

PRI4 

(R53 1 -R670)/(R53 1 -R670) 

0.85; p<0.01 

0.63; p<0.05 

L 

NDW1 

(R870-R 1 240)/(R870+R 1 240) 

0.76; p<0.05 

0.60, ns 

L 

NDVIuindwi 

(NIR-RED)f(NIR + RED ** 

0.19 : ns 

048. ns 

L 

b. Hardwoods (HW). Duke. NC (n=7) 

PRI4 

( R53 l-R670V(R53 1 -R670) 

0.95; p<0.0l 

0.48, ns 

NL 

Dmax 

Dmax= Highest 1 st Derivative, 650-750 nm 

0.83; p<G.0L 

0.40, ns 

NL 

ND1I 

(R820-R1 650 V ( RS20+R 1 650) 

0.79; p<0.01 

0 .34, ns 

L 

EVIbamlm 

2.5(NIR-RED}/(NIR+6RED-7.5BL UE+ if 

0,84 p<0.01 

0.41, ns 

L 

ND l / Uimtsat 

(NIR-Red)f(NIR+RED) 

0.63 : p<0.05 

0.19, ns 

L 

c. Tall Grassland fTG). Konza. KS (n=14 L eombinina 1 & 4 year bum} 

G94 

(R30Q-R70Q)/{ R800+R700) 

0.95; pcQ.OI 

0.90; pO.QI 

NL 

Dmax 

Maximum of the I &1 Derivative in 650-750 nm 

0.93; d<0.0 1 

0.95; p<0 0l 

L 

NDW1 

( RS70- R 1 240 V( RS70+R 1 240) 

0.86; p<0 01 

0.9I; p<0.01 

NL 

Ei'f Landsat 

2J(RNIR~Rred)/(RNW+ 6Ried- 7 5RbIue+ if 

0.83; p<0.01 

0.89; p<0.01 

NL 

ND ViLandwt 

(NIR-Red)/(NIR+RED) 

0.63 r ns 

0.67; p<0.05 

L 

d. Miomho Wtiodland (MW), Monau. Zambia (n=24} 

G32 

R750, 700, 450 

0.83; p<0.0l 

0.8I; p<0.01 

NL 

Dmax 

Maximum of the I st Derivative in 650-750 urn 

0.77; p<0.0l 

0.87; p<0.01 

NL 

Dmax/D704 

Dmax/D704 

0.79; p<0.Q I 

0.80; p<0.01 

NL 

Landsat 

2. 5(RNIR-Rred)/(RNlR+6Rred- 7.5Rbitte+ 1) 

0.73; p<0.0 

0.88; p<0. 01 

L 

ND l 

(NM-Red)7(NIR+RED) 

0.52 , ns 

0,60; p<0.05 

NL 


*NEP (net ecosystem production), GEP (gross ecosystem production), L (linear relation), NL (non-linear regression) 
**NIR (average R760-900 nm), Red (average R620-690 nm), Blue (average R450-500 nm), after the Landsat bands 


results illustrate the possibility for establishing a common 
spectral approach based on high resolution Vis for tracking 
C0 2 ecosystem dynamics. Dmax was able to provide separa- 
tion for most of the data, but it saturated at Dmax > 0.65 for 
NEP > 25 //mol m -2 s -1 (Fig. 8). The use of derivatives, 
reducing the effects of variable irradiance, canopy structure 
and bidirectional properties, offers an important potential when 
analyzing reflectance time series. 

It is important to note that the Landsat broadband Vis never 
performed as well as the narrow-band indices. At the HW, TG 
and MW sites EVI performed almost as well as the “optimal” 
narrow-band Vis, but not at the LP site (Table VI) or when an- 
alyzing the combined data set (Table VII). 

IV. Discussion 

We were able to evaluate Hyperion reflectance stability over 
the full instrument spectrum and large portion of the EO-1 life- 


time, for which there was available data. Since the collection of 
sufficient number of free of clouds images to generate the time 
series took multiple years, we constructed a seasonal time line 
by arranging the data consecutively by acquisition day from any 
given year. 

Our results from the three PICS demonstrate the spectral sta- 
bility of Hyperion’s reflectance time series, and provide spa- 
tial uniformity descriptors at each of these well-known sites 
where the time series have been collected. This spectral stability 
verifies the potential of Hyperion’s reflectance time series for 
ecosystem monitoring. As new images are added to the seasonal 
time lines, the ability of the reflectance time series to assess veg- 
etation function will likely improve. Future work is required to 
refine the atmospheric correction process and to improve the re- 
flectance estimates in spectral regions with lower stability (e.g. 
2000-2450 nm region, and near water absorption bands). Fur- 
ther improvements in the atmospheric correction approach, the 
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TABLE VII 

Spectral Bio-Indicators of Vegetation Function With Best Correlation to C0 2 Uptake for Multiple Sites (Combined Data Set From 

the Three VFS, n = 52) 


Spectral 

indicator 

Formula 

r 2 toNEP + (L/NL)* 

iNoGEPJfL/NL) 

Dntax 

Highest 1 st Derivative, 650-750 nm 

0.73. p<0.0! (L) 

0.77 . dO.OI (L) 

DP22 

Dmax/D (ina![+12) 

0.65, p<0.0l (L) 

0.74. p<0.01 (NL) 

NDW1 

R(37G-I240)/R( 870- 1240) 

0.74. p<0.01 (NL) 

0.67. pcO.Ol, (NL) 

MCARIa 

|(R700-R675)-(0.2*(R700-R550))*R(700/R675)| 

0.41,p<0.01 (L) 

0,75, p<Q.Gi (L) 

PR14 

(R53 l-R670)/(R53 1 -R670) 

0.66, p<0.01 (NL) 

0.62, p<0.05 (NL) 

ADI I faitjdsaf 

(NIR-R)/(NIR+R)** 

0.56, ns (NL) 

0.59. ns (NL) 


t NEP (net ecosystem production) 

* GEP (gross ecosystem production) 
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Fig. 6. Relationship of Hyperion’s best performing Vis in relation to NEP at each study site (p < 0.01). At Duke (Fig. 6(a)), the photochemical reflectance 
index PRI4 using green and red wavelengths (PRI4 = (R531 — R670)/(R531 — R670), performed well at both LP (r 2 = 0.85) and HW (r 2 = 0.94) 
forests. However, the index saturated at PRI4 ~ 0.2, not being able to effectively discriminate NEP > 20 //mol m -2 s _1 . At the TG Konza site, the most 
successful VI was narrowband G94 (r 2 = 0.95, Fig. 6(b); G94 = (R800 — R700) / (R800 + R700)), which saturated at G94 ~ 0.7, not being able to separate 
NEP > 27 //mol m~ 2 s _1 . At the Mongu MW site, the most successful VI was G32 (r 2 = 0.83, Fig. 6(c), G32 = (R750 — R445)/(R700 — R445)) combining 
wavelengths from the blue, red and NIR regions. At low values of G32 < 0.3 this parameter saturated, not allowing separation of NEP < 5 //mol m -2 s _1 . 


incorporation of topographic information, and additional cali- 
bration of Hyperion reflectances to targeted in situ field spec- 
tral measurements, are required to expand the number of useful 
reflectance bands at VFS. These improvements are especially 
needed in the 850 to 2500 nm range where they may allow better 
reflectance retrieval and thus, narrow-band Vis using bands in 
these ranges may be better correlated to NEP and GEP than in 
the current analysis. 

Hyperion’s long-term reflectance measurements offer the 
potential for cross-comparison with other optical sensors. 
Cross-calibration of current multispectral sensors to hy- 
perspectral instruments such as Hyperion, could serve for 
inter-comparisons and retroactive data processing. Currently 
there is a consistent NDVI record that extends more than 2 
decades, generated from the Advanced Very High Resolution 
Radiometer (AVHRR) instruments, which provides a histor- 
ical perspective on vegetation dynamics necessary for global 
change research. Long term spectral records are also avail- 
able from SPOT- Vegetation, SeaWiFS, Moderate Resolution 
Imaging Spectroradiometer, Landsat ETM+ data and others. 


While the spectral vegetation indices (Vis) are considered more 
directly transferable between sensors, and an intercalibrated, 
sensor-independent NDVI data record could be feasible, there 
are significant challenges and higher discrepancies associated 
with the inter-comparison of biophysical products, which are 
generated using different instruments with varying spectral 
and spatial resolution, modeling approaches and ground data 
[50] — [52]. 

Simultaneous calibration of old and new sensors against each 
other over a common set of targets on the Earth is a key re- 
quirement for maintaining long-term records (http://wgcv.ceos. 
org/). Our results indicate that the use of space-borne spectrom- 
eters allow assessment of the spectral stability and spatial vari- 
ability of the calibration sites. The reflectance measurements at 
the PICS confirm that they meet the CEOS/WGCV surface re- 
flectance requirement for above 30% reflectivity and flat spec- 
tral signature over the 600-2000 nm range. 

In recent years, a growing number of studies have recognized 
the importance of collecting spectroscopy observations at flux 
towers to better understand the relationship between vegeta- 
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Fig. 7. Temporal variation in the best performing spectral indices for Konza, Duke (LP), Mongu (reflectance means (x) by DOY). Dotted line is NEP while dashed 
line is precipitation. 
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Fig. 8. Relationship of Hyperion’s derivative index Dmax (the maximum 
amount of the first derivative in the 650-750 nm region) to the C0 2 as- 
similation (NEP, //mol m -2 s _1 ) across the three VFS. Dmax had closest 
relationship when using the combined data set (n = 52, p < 0.01). However, 
this parameter saturated at Dmax > 0.65, precluding the separation of 
NEP > 20 //mol m -2 s -1 . 


tion function and optical signals [53]. Our results at the three 
VFSs demonstrate the close link between Hyperion’s Vis and 
carbon flux measurements and the potential for scaling the tower 
measurements to local, regional and potentially global levels. 
Multiple observations at individual sites are required to test the 
ability of the Vis to describe seasonal and transient variations 
in fluxes. It is noteworthy that the Vis with stronger relation- 
ships to the C0 2 parameters were derived using continuous re- 
flectance spectra or reflectance derivatives, and included wave- 
lengths associated with chlorophyll content and/or chlorophyll 
fluorescence (e.g. DPI, Dmax/D704 and PRi). Since these in- 
dices cannot be calculated from broadband multispectral instru- 
ment data, the opportunity to exploit these spectrometer-based 
Vis in the future will depend on the launch of satellites such as 
En MAP and HyspIRI. 

Space borne spectrometers offer potential for establishing a 
generalized approach for scaling C0 2 dynamics across multiple 


sites and to local and regional levels. Site-specific equations 
provided more accurate results than the multisite groupings, 
although the trade-offs between the improved accuracy of a 
site-specific algorithm versus an algorithm for general use have 
yet been established. The development of common (i.e. global) 
spectral approaches to track vegetation function and scale up 
carbon estimates requires diverse spectral data sets, comprised 
of the major ecosystem types in a broad range in functional 
states. 

The performance of narrow-band Vis from satellite time se- 
ries remains to be thoroughly examined for various vegetation 
types exhibiting a range of local conditions. Future work will ex- 
pand the analysis by extending the time series at our study sites 
as well as adding additional sites representing other ecosystem 
types, including rain forest, temperate and sub-arctic vegeta- 
tion. Testing algorithms over longer time series and more di- 
verse cases will increase our ability to determine fluxes with 
higher confidence. Using hyperspectral data operationally will 
allow us to develop new more dynamic and accurate approaches 
to predict ecosystem parameters. 

V. Conclusions 

With over 1 1 years of data record including more than 65,000 
globally distributed hyperspectral images collected to date (as of 
July 20 2012), the Hyperion data record provides an important 
and unique resource for various Earth observations. This study 
provides a step towards better understanding and broader use 
of the Hyperion reflectance time series for monitoring of vege- 
tation function and ecosystem dynamics. Hyperion reflectance 
measurements were shown to be stable over time in all spec- 
tral regions. The most stable portion of the spectrum was the 
visible and near infrared region, providing higher confidence 
for common use of chlorophyll based Vis. In three different 
ecosystems, Hyperion’s Vis were found to be correlated with 
CO 2 flux parameters. The Vis with the strongest relationships 
required continuous spectra, or numerous wavelengths associ- 
ated with chlorophyll content and/or first derivative parameters, 
which are not available from broad-band multispectral instru- 
ments. This study showed their use with Hyperion’s time series, 
demonstrating their close relationship to ecosystem carbon dy- 
namics and potential for improved monitoring of carbon fluxes. 
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